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Abstract 

If dark matter (DM) has non-zero direct or transition, electric or magnetic dipole 
moment then it can scatter nucleons electromagnetically in direct detection experi- 
ments. Using the results from experiments like XENON, CDMS, DAMA and CO- 
GENT we put bounds on the electric and magnetic dipole moments of DM. If DM 
consists of Dirac fermions with direct dipole moments, then DM of mass less than 10 
GeV is consistent with the DAMA signal and with null results of other experiments. 
If on the other hand DM consists of Majorana fermions then they can have only non- 
zero transition moments between different mass eigenstates. We find that Majorana 
fermions with mass > 38 GeV and mass splitting of the order of (50-200) keV 
can explain the DAMA signal and the null observations from other experiments and in 
addition give the observed relic density of DM by dipole-mediated annihilation. This 
parameter space for the mass and for dipole moments is allowed by limits from L3 but 
may have observable signals at LHC. 



1 Introduction 



Experimental observations mainly of dynamics of spiral galaxies and galaxy clusters indicate 
the existence of dark matter (DM). Cosmological observations confirm the existence of DM 
and in addition show that the bulk of it must be non-baryonic py. In this paper we consider 
weakly interacting massive particles (WIMPs) as candidates for DM, but we adopt a model 
independent and phenomenological approach. 

The mass and cross section of the DM (which is expected to have a local density of about 
0.3 GeV/cm^ and velocity w.r.t the Earth of about 200 km/sec [2]) is probed by direct 
detection experiments like XENON [3J, CDMS [1], DAMA [5j and COGENT |6j. These 
experiments detect DM scattering off nuclei by measuring the recoil energy of the nuclei. 
The energy threshold of such detectors is typically of the order of a few keV. Of these 
experiments, the DAMA experiment observed an annual modulation in its signal which 
could have been due to DM scattering. However none of the other experiments conducting 
direct DM searches have seen evidence of such an event. Theoretically if one considers spin- 
independent interaction [3, |H] of these WIMPs with the nuclei then it is found that only in 
the low mass range from about 5 to 10 GeV one can reconcile DAMA signal with the null 
results from other experiments. There have been proposals of inelastic scattering of DM [9] 
by nuclei which results in heavier DM mass being allowed for explaining all the existing data. 

Although DM has zero electric charge it may couple to photons through loops in the form 
of electric and magnetic dipole moments. Here we study WIMPs which are endowed with 
such dipole moments and thus can interact feebly via electromagnetic interaction [TUt [TT]. 

We first will consider the case that DM is a Dirac fermion. The effective Lagrangian for 
coupling of a Dirac fermion x having an electric dipole moment V and a magnetic dipole 
moment /i to a electromagnetic field J-''^'^ is 

Celast = wAt^ + 1^V)XJ''"'. (1) 

In this case, we have elastic scattering of DM by nuclei through a photon exchange 
(FigH]). This is studied in detail for the case of electric dipole moment in Section [2] and of 
magnetic dipole moment in Section HI 

Next we consider the case of DM being Majorana fermions. These have only non-zero 
transition moments between different mass eigenstates. Their interaction with photons is 
described by 

C-inel = -- X2Cri^u{fl'12 + 75^^12) Xl (2) 

where /ii2 is the transition magnetic moment and T>i2 is the transition electric moment. 
Majorana DM would have inelastic scattering off nuclei through a photon exchange. We 
study the case of transition electric dipole moment in Section [3] and of transition magnetic 
dipole moment in Section HI 
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Figure 1: Electromagnetic scattering of a proton with DM with non-zero dipole moments. 



In Section Owe calculate the bounds on the dipole moments which would give the desired 
relic density abundance of DM. 

In this paper we present two main types of results: 

1. We find bounds on dipolar moments, direct and transition, coming from WIMP search 
experiments, which update the results in [TT]. In addition we show that there are 
bounds on dipole moments from single photon search at LEP, an effect not discussed 
in [H]. 

2. We find regions in parameter space that are consistent with the positive signal from 
DAMA and with the null results from other experiments. 

All these results are presented in Section [61 In Appendix A we give some details of our 
DM scattering calculations, and in Appendix B we show the DM annihilation cross section 
calculations which are relevant for the relic density results. In Appendix C we review the 
contact scalar interaction case, which is useful for comparison with our results. 



2 Electric Dipole Moment Interaction of Dark Matter 

The differential cross section for a DM-proton elastic scattering via interaction of the electric 
dipole moment of DM with the proton charge is given by 



dEn Attv^Er 

where V is defined in ([T]). Here Er is the recoil energy of proton and v is the speed of DM 
relative to the nucleus. The differential rate for nuclear scattering of DM through electric 
dipole moment interaction for a nucleus of Z protons and DM of mass is given by 

'""-Z^Nr^f mv4^dv (4) 



dER Jv>v^,„ dEn 
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Figure 2: Plot shows the allowed regions for DM electric dipole moment with varying DM 
mass for elastic scattering for different experiments. Shaded region shows the allowed pa- 
rameter space for DAMA which is consistent with all other experiments. 

where iVy is the number of target nuclei in the detector and is the local DM density. 
The velocity distribution f{v) and minimum speed of DM fmm, for a given energy threshold 
ERmin are given by p| 



where rriM is the mass of the target nucleus, /itv = 'm^rriM / {m^ + is the reduced mass of 
the DM-nucleus system and Vq = 220 km/sec. The factor appearing in eqn.(jl]) shows the 
fact that this is a coherent scattering of the Z protons in the target nucleus. To determine 
the expected event rate R we integrate eqn. (jll) over the nuclear recoil energy: 




(5) 



(6) 




(7) 



Here e{QEji) is the efficiency of the detector which depends on the recoil energy and 
Q is the quenching factor that relates the observed energy with the actual recoil energy 



Experiment Target nucleus Quenching 



factor((5) 


DAMA 


Na 


0.3 




I 


0.09 


CDMS 


Ge 


1 


XENON 


Xe 


1 


COGENT 


Ge 


0.2 



Table 1: Quenching factors for different experiments used in this analysis. 

i.e Eobs = QEr. Q depends on the target nucleus and the nature of the detector. Ei 
and E2 are limits of the observed energy interval in a direct search experiment. Therefore 
Ei/Q and E2/Q give the corresponding limits for the recoil energy Er. Observed energies 
are usually quoted in units of electron equivalent energies (keVee). The detectors usually 
detect scintillation caused by electrons and thereby measure energy transferred to electrons 
by the nucleus. Q then represents the efficiency with which the recoil energy of nucleus is 
transferred to electrons which are detected by scintillation. The quenching factors for the 
different experiments analyzed here are given in Table [H 

In Figj2] we show our results. The experiments give an upper limit which is a function 
of the DM mass m^, except for DAMA whose result is represented as an allowed band. We 
would like to point out that one can get a bound from the e+e' XXI process whose signal 
is a single photon detection with missing energy. From the analysis of the collaboration L3 
p!2] we get V < 6.6 x 10~^^ e-cm, valid for < 38 GeV; this limit is also shown in FigEl 
We shall see the implications of this Figure in Section El there we shall also discuss the issue 
of the DM relic density. 

3 Inelastic Dark Matter and Electric Dipole Moment 
Interaction 

In this Section we consider the case that the WIMP is a Majorana particle, with a a transition 
electric dipole moment P12 as defined in ([2]). In this case we can have inelastic scattering 
Xi + N ^ X2 + N where xi ^-^id X2 are two different mass eigenstates, and in general there 
is a mass difference between xi arid X2, S = 1112 — nii. Due to this mass difference, the 
minimum DM kinetic energy needed for the nucleon scattering becomes higher. DAMA has 
lower detection threshold compared to most experiments and therefore can be more sensitive 
to this scattering mode than the other experiments. This was first proposed in 

If we drop terms of higher order in 5 (compared to we end up with eqn.Q for the 
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Figure 3: Plot of DM electric dipole moment against the mass difference 6 for inelastic 
scattering for different experiments. Shaded region shows the allowed parameter space for 
DAMA which is consistent with all other experiments and in case of > 38 GeV it is also 
consistent with the relic density. 
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Figure 4: Plot of allowed DM mass and mass splitting 6 for a fixed V which gives the correct 
relic density of DM. 



differential cross section. That is, 

d(T inelastic d(Jelastic 

dEn dEji 

However the total cross section in the elastic and inelastic case are related to each other 
by a factor which enforces the condition of the minimum kinetic energy required for inelastic 
scattering. This relation has been derived at the end of Appendix A, and is given by, 



'^inelastic ^^elastic- (9) 

But we do not use the total cross section in our analysis, only the differential cross 
section and therefore the square root factor does not figure in our calculations. We integrate 
the differential cross section per unit recoil energy over the energy interval in which the 
experimental data have been observed. The only change from the analysis in the elastic case 
is that the minimum speed for scattering is now given by 



+ b\ (10) 
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where notation has the same meaning as in the previous Section. 

In FigO] we plot our results. We show a line for each direct search experiment and the 
allowed range is on the right of the line, except for DAMA that again is an allowed band. 
In the Figure we also plot the limit coming from LEP mentioned before, 1^12 < 6.6 x 10~^^ 
e-cm, which is the same than in the direct transition case if \6\ <C m^. The collider limit 
applies provided < 38 GeV. 

Since in the inelastic case we have three parameters, namely, X>i2,m^, and 6, it is in- 
structive to present FigH] where we fix and show the constraint in the {m^, 6) plane. 
Here the allowed regions is on the right of the curves of the different experiments except for 
DAMA which is again an allowed band. 

We shall discuss the implications of all these limits in Section [6] and especially for the 
relic density in Section O 

4 Magnetic Dipole Moment Interaction 

The differential cross section per unit energy transfer for elastic scattering by a magnetic 
dipole moment interaction is given by 



In the case of inelastic scattering we get the same formula with /x substituted by the 
transition magnetic moment /ii2, provided we drop terms of higher order in the mass differ- 
ence 6. Again, the only difference between elastic and inelastic is in the kinematics, because 
the minimum speed for inelastic scattering is given by eqn. lJTOl) . Fig I5] and Figj6] show the 
results for magnetic dipole moment which are essentially the same as those for electric dipole 
moment in terms of the allowed parameter space from various experiments. In FigJT] we plot 
the allowed regions in the space {m^, fiu)- 

In the next section we discuss the implications of dipolar DM vis a vis the DM relic 
density. 

5 Relic Density of Dipolar DM 

The relic abundance of dipolar DM is determined by the annihilation cross section xx ~^ ff 
through a dipole vertex. In the case of Majorana fermions only non-identical fermions can 
annihilate {xiX2 ff) through the dipole channel as only the transition dipole moments 
(Pi2,/ii2) are non-zero. When the mass difference 6 between xi and X2 is small the cross 
section for Majorana annihilation process is identical to that of the Dirac fermions (with 
V, n replaced with V12, 1^12) ■ 




(11) 
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Figure 5: Plot shows the allowed regions for DM magnetic dipole moment with varying 
DM mass for elastic scattering for different experiments. Shaded region shows the allowed 
parameter space for DAMA which is consistent with all other experiments. 

The annihilation cross section for electric dipole annihilation of DM for xx ~^ ff is given 
by (see appendix B) 



where v = 2y 1 ^ is the relative velocity of the two annihilating WIMPs. The thermal 

averaged cross section, (aVrei), can be parameterized in terms of the temperature, T = 
m^(f^)/3, where v = Ure//2, such that 




(12) 



and the expression for the magnetic dipole annhilation case is 




(13) 





(14) 
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Figure 6: Plot of DM magnetic dipole moment against the mass difference S for inelastic 
scattering for different experiments. Shaded region shows the allowed parameter space for 
DAMA which is consistent with all other experiments and in case of > 38 GeV it is also 
consistent with the relic density. 
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Figure 7: Plot of allowed DM mass and mass splitting 5 for a fixed /x which gives the correct 
relic density of DM. 



where n and the parameter ao for magnetic and electric dipole interactions take the following 
values 



Magnetic Dipole: ctq 
Electric Dipole: ctq 



eV 
47r ' 

487r 



n = 



n 



(15) 
(16) 



Considering cold dark matter (CDM) we want to estimate the relic density in the context 



of dipolar interactions of such CDM candidates. The expression for Xf 



, where the 



/ 



subscript / indicates the freeze out condition, is given up to a reasonable approximation by 



- I n + i j ln{ln[0.038(n + l){g / gll^)mpim^a^]} 



;i7) 



where mpi is Planck mass, while g and g^^s are the effective number of relativistic particles 
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at the time of decoupling. The DM rehc density is then given by [T^ 



nh^ = 0.34 ^ . 18 



The hmit on cold dark matter density from WMAP [13] is fi^n/i^ = 0.1099 ± 0.0062. 
Using ( JTSll . this gives us the bound on T>, which gives the acceptable relic density of dark 
matter, to be ~ 2.5 x 10"^^ e-cm. For the magnetic dipole case again using (ITSll we find that 
the dipole moment must be ~ 8.2 x lO-^/i^- 

For the case of non-identical Majorana annihilation the annihilation cross sections are 
identical to the case of Dirac fermions and the limits on the Dirac electric and magnetic 
moments from relic density abundance are identical to the bounds on transition magnetic 
moments given above. 



6 Results and Conclusions 
6.1 Direct dipole moments 

We have explored the parameter space of electric and magnetic dipole moments of DM from 
the results of nuclear recoil experiments. We find that the limit on direct electric dipole 
moment is P < 1.6 x 10"^^ e-cm if the DM mass > 10 GeV, see Figj2l In this mass 
range there is no allowed parameter space which is consistent with the positive signal from 
DAMA and with the null results of other experiments. 

For a DM mass = (3 — 13) GeV and a moment V ~ 2.5 x 10"^° e-cm there is a region 
in parameter space where the DAMA signal is not ruled out by other experiments, see Figj2l 

For smaller masses, < 3 GeV, the Figure shows that COGENT rules out DAMA. In 
Appendix C we review the scalar interaction mediated scattering (see FiglH]) where we see 
that DAMA is allowed by COGENT for any DM mass below about 15 GeV. This can be 
attributed to the fact that the DM interaction via electric dipole moment has a factor of 
recoil energy [En in the denominator of its cross section in eqn.(jl]). This means that for low 
threshold this cross section is large and in the case of COGENT this threshold happens to 
be smaller than that for DAMA. But as the DM mass increases this advantage in favour of 
COGENT is nullified by the DM mass factor which again appears in the denominator of the 
cross section formulae. 

The limit on direct magnetic dipole moment is /i < 10~^^b if the DM mass > 25 
GeV, see FigJSl In this mass range there is no allowed parameter space which is consistent 
with a positive signal from DAMA with null results of other experiments. For DM masses 
rriy. = (3 — 12) GeV and /i ~ 6.3 x lO^'^/xs there is a region in parameter space where the 
DAMA signal is not ruled out by other experiments, see FigO Actually, in the low DM 
mass range the collider bound [12j is more restrictive than the direct WIMP search limits, 
as can be see in FigJSl 
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If the only interaction of DM with standard model particles had been via electromagnetic 
dipole interaction then the relic density of the DM particles would be determined uniquely 
from the annihilation cross section of the EM-dipole mediated process XX ~^ ff ^ind m^. 
Such a calculation shows [H] that to get relic cold DM density consistent with WMAP 
measurement the dipole moments have to be in the range of P ~ 5 x 10"^'' e-cm or /i ~ 
10~^fiB when = (10 — 1000) GeV. The limits on DM dipole moments that we obtain 
on the basis of nuclear scattering experiments are much lower. This means that the DM 
particles must have some other more dominant interaction which decouples at a temperature 
Tf ~ m^/10 to give the correct relic abundance while the dipole interactions decouple much 
earlier and have no bearing on the present DM density. These other dominant interactions 
however must not dominate over dipole interactions in nuclear scattering in order for our 
bounds on dipole moments to be meaningful. A good example of such a situation is seen in 
certain classes of WIMP models [2], where the interaction with standard model particles can 
be by Higgs exchange which gives rise to spin-independent (SI) interactions with nucleons or 
by Z exchange which gives rise to spin-dependent (SD) interactions [11]. When the nuclear 
scattering is by SD interaction then there is no coherent enhancement of the cross section 
by the atomic number. The bounds on the cross section for SI interactions is therefore more 
stringent than the bounds on SD interactions. The bounds for SI interactions from nuclear 
scattering experiments is asi < 5 x 10~^^ cm^ and for SD interactions it is ctsd < 5 x 10~^^ 
cm^ (for ~ 100 GeV). Our bounds are relevant for WIMPs which can have a large spin- 
dependent cross section aso — 0.3 x 10~^^ cm^ to give the correct relic density, but this cross 
section is too small to be observed in nuclear scattering experiments. 

6.2 Transition dipole moments 

In FigsJS] and [6] we plot the constraints on the transition dipole moments V12 and fiu, 
respectively, versus the DM mass m^. We have fixed three values for the mass = 14,40 
and 100 GeV, in order to illustrate our findings. 

We see that a common feature of the plots is that the region where DAMA is consistent 
with other experiments is quite large and extends to higher and higher inelasticities. However 
for < 38 GeV the region is further constrained by the collider bound from L3, and the 
allowed region is reduced. The allowed area is more and more reduced as we lower and 
for < 10 GeV there is no region at all. For > 38 GeV there is no collider bound and 
the region is large. 

The most interesting aspect of the experimentally allowed parameter space for transi- 
tion dipole moments is that the cosmologically preferred value (for getting the correct relic 
abundance) Vij ~ 10~^^ e-cm or fiij ~ lO'^fts [llj is consistent with the allowed values 
DAMA and other experiments when > 38 GeV and the mass split 6 ~ (50 — 100) keV. 
Thus, we find that transition dipole scattering has a large parameter space where the results 
of DAMA, the null results from other experiments and cosmological relic abundance are all 
consistent. 
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The event rates depend upon and P or as V^/ rriy^ or /x^/ respectively. The dipole 
moments can increase by four orders of magnitude from the V — 10"^° e-cm to P = 10"^^ 
e-cm (and similarly for the case of magnetic moments) and the same event rates would be 
obtained if mass is increased by eight orders of magnitude from 10^ GeV to 10^'^ GcV . 
The range of DM mass for which DAMA is consistent with cosmology and L3 is therefore 
quite large, 38 GeV < < 10^° Gev. This parameter space of and electric /magnetic 
moments can give rise to a missing energy signal via the process pp — > XiX27+hadronic jets, 
and the dipolar model of dark matter may be testable at the LHC. 
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Appendix A: Cross Section Calculation for Electric and 
Magnetic Dipole Moment Interactions with Nuclei 

Here we consider the elastic scattering process X + p — * X + We calculate the cross 
section in the lab frame where the proton is at rest initially. We assume that even after the 
scattering the proton remains at rest since the momentum transferred is very small in this 
process. The entire calculation is done in the non-relativistic limit. 

Electric Dipole Moment 

The initial momenta of DM and proton are denoted by ki and Pi respectively, while the 
final momenta are denoted hy kf and pf respectively. The amplitude squared for this process 
is given by 

= ^[64(g • Pi){h ■ Pi) + 32mJ(g • pi) - 16(g • pif - %q\q ■ pi) - 64:{ki ■ pif] (Al) 

where q is the momentum transferred, while and rUp are masses of DM and proton 
respectively. In writing the above equation we have made use oi kf — ki — q and pf = Pi + q. 
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Here, (g ■ Pi) = vipEji and ki ■ Pi ^ m^nip. Hence we can drop terms containing (g ■ pi) and 
its higher orders, keeping just the last term in eqn. OAip . This gives us 



16e2p2 



\M\' = -^Km,)^ (A2) 

Using = E"^ — 2mpE^ ^ —2mpEji we finally arrive at the following expression for the 
amplitude squared. 



8e'V 



|2 — rrj^ 

Now, the differential cross section is given by 



|A<P = -^™>p- (A3) 



- iSk. <'^>* ''<^- + - - i-^i^ 

1 

L5{mp + E,-Ep-E^)\M\^ 



GATT'^mpEiv E2E. 
1 Ifc 



(^dQ\M\^ (A4) 



QAn'^mpEiV Ej 

where Ei and £^2 are energies of the initial and final DM states. Also the recoil energy can 
be written as follows: 

1^^ \ki\'^ + \kf\^ -2\ki\\kf\ cos 6 

Er = = 

2mp 2mp 



therefore, 

dER = - '•"'"•"■^' rf(cosg). 
nip 

And we know that dQ = —2'n'd{cos6), so we can write 



\ki\\kf\ 



dn = 27r S dER, (A5) 



Using this in eqn. (lA4p we get 

da 



dER 327tE,v Ep\k,\ 



327iEfEpV^ 

We use the approximations, Ef ^ and Ep ^ nip, which gives us 

dcr 1 



\M\^. (A6) 



dER 327rm?,m„f2 



\M\^. (A7) 
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Substituting for from eqn. flASp we get 

da e^V 



which is same as eqn.(l3]). 



Magnetic Dipole Moment 



(A8) 



The trace in this case is same as in the case of EDM but without any 75 in it. The 
amphtude squared for this process is given by 

= — — [4(g ■ Pi){ki ■ Pi) - 4(g ■ Pi)m\ - (q ■ pif - A{ki ■ pif + Am^ml] 



[4{q-piYml] . 



Here we drop terms which are higher order in g ■ pj except the last term which is enhanced 
by a factor of in the denominator. Thus we get 



\M\^ 
where 



q^ 



16eV [/ \ ( \ 2 (u 2 2 {q-Vi?ml 

[q ■ Pi){ki ■ Pi) - {q ■ Pi)m^- [ki ■ Pi) + m^nip - ■ 



(AlO) 



q-Pi = rUpER 
ki-pi^ m^rrip 
{ki-Pif 



2 2 I 222 



Using these results in eqn. flAlOp we get 



16e_2/i' 
q' 



m^rripER - mpm^^En - m^mlv'^ - 



Again using q^ ^ —2mpEji gives 



8eV 
8eV 

Er 



m^mpV^ + Er — m^mpER 



mlmlEl 



m\ER 



ttiIEr 



m1m„v'^ H — 



m^nipER 



(All) 



(A12) 
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The differential cross section per unit recoil energy is given by eqn. flATp . Substituting 
from the above equation in eqn. flA7p we get 



da 




1 + 


Er 


Er 1 


dEn 




2mpf 2 








1 + 


Er 


^Er 1 




4=7tEr 




2m^f 2 



where fip 



my/rrip 
+ nip 



is the reduced mass of the DM-proton system. The last term in the 



above equation can be dropped for large values of DM mass. And so we get the following 
expression for the differential cross section per unit energy, 



da 



1 + 



E 



R 



dER 4:71: Er [ 2/ipt;2 
Total Cross Section for Inelastic Scattering 



(A13) 



In the case of inelastic scattering we consider the WIMPs as Majorana particles. As 
mentioned before in the main section that the differential scattering cross section in both 
elastic and inelastic scattering is the same. But the total cross sections in the two cases are 
related to each other as in eqn.®. Now for inelastic scattering we have, 

kl — kl = 2yUjv5 
or kl = kl — 2fiiyS. 

The magnitude of the momentum transferred (working in CM frame) is then, 

- -1 ^ -2 
2kl ( 1 



- kf + kl- 2k'),k'), cos 6 



26 



cos 9 \ — 2fiN5 



where we have used kf = nj^v . And so 



dq' 



But dq"^ = 2m N dER, therefore 




26 



■(i(cosi 



dE 



R 



kl 
TUN 



26 



■(i(cosi 
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The above equation in the elastic case becomes 

dEn = ^-dicosO). 

rriN 

Now the total cross section for elastic scattering is given by, 

/d<7 elastic , 7-1 
~dE — 



27rkl f da elastic 



I 



mjv J dO, 

For the inelastic case we have, 



(i(cos 9). 



2'Kkl 2S f daelasticj, ri\ 
(^inelastic = Wl ^ / 17^ d(cOS 6/) 

fTlN V I^^nv J d\l 



-L 9 (^elastic ■ 

Appendix B: Annihilation Cross Section of Majorana 
DM 

The process involved here is Xi(^i)X2(^2) f{Pi)f{P2), f being a fermion of charge e with 
negligible mass. Now we work in the CM frame where, 

ki^{Ei,0,0,k), k2^{E2,0,0,-k), 
Pi — {Ef, Ef sin 9, 0, Ef cos 9), pi — {Ef, Ef sin 9, 0, Ef cos 9), 

? = Pi+P2, q^^AEj. 

The amplitude squared for this process is given by 



2 2 

l-Mr = ^g"/b(A;2)t'(A;2)r^aM(A;i)M(A;i)r,;3][w(pi)w(j>i)7^t'(P2)t^(j>2)7l (Bl) 

where Ff^ce — CT/^a and g — /i for magnetic dipole interaction, whereas for electric dipole 
interaction — '^ixalb a-nd g — T>. For Majorana particles. 
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In the Majorana representation of 7-matrices, C = —70 and the Majorana condition 
imphes that ip = ip'^ = Cip'^ = ip*. Imposing ijj = ijj* on (IB2p we get the relation Ug = v*. 
Using this relation we can rewrite eqn. llBip as 

2 2 

l-^l' = ^g"/b(fc2)u^(A;2)7or^„n(A;i)t;^(fci)7or,^][n(pi)n(pi)7^t;(p2)t)(p2)7l- (B3) 



Now the spin sum rules for Majorana fermions are given by [16] 

J2 ^s{p)vj{p) = (/+ m)C^ (B4) 

s 

J2^sip)u^{p) = i^-m)C^ (B5) 

s 

where = 70. We then take the sum over spins of eqn. (lB3p and use eqns. (lB4p and (IBSp 
to obtain 

2 2 

= ^Q''Q^Tr[i^2 - m2)r^„(^i + mi)r,^] Tr[^i7^;^27l (B6) 

where m\ and m2 are masses of Xi and X2 respectively. Now we discuss the two cases, 
beginning with the electric dipole case first. 

Electric Dipole interaction 

Eqn.( 1B6p with g = V and r^,^ = (T^i/75 yields the following expression for the averaged 
amplitude squared, 

4g4 



= —r^{Pi ■ P2)[{h- Pi){ki ■P2) + {h ■ Pi){k2 ■ Pi) - {Pi ■P2)mim2] (B7) 



where q = Pi + P2- Now in CM frame, 

ki ■ pi = EiEf — kEf cos 6 ^1 ■ P2 = EiEj + /ci?^ cos 6 

k2 ■ pi = E2Ef + kEf cos 9 k2 ■ P2 = -E2-E/ — kEf cos 6' 

Pi • P2 = Y q =AEf 

Using these in eqn.( ]B7p we get 

JM^ = Ae^V'^[EiE2 - k^ cos^ 9 - rxixm-^ 

= Ae^V^E,E2 ^) + (l- i^-^'^-^D] ,os^0 

\_\ E1E2 J \ 2E1E2 

s ^ ^ (m? + rrin) , ,^ ^ 
where we use A; = - - £^1^2 - ^ — ^ with s = {Ei + £2) • 



:B8) 



The differential cross section in the CM frame is given by 
da _ 1 Ef 



dQ 64:n'^EiE2Vrei \/s 



EIE2 J 



[s — m\ — m^) 
2E1E2 



cos^e 



(B9) 



where Ef = Integrating above equation over dD. gives the total cross section as 



aVrel 



in 



^ mim2\ 1 A {s — mf — 7712) 
E1E2 J 3 \ 2E1E2 



(BIO) 



Eqn.( |T2l) can be obtained by taking mi = m2 and Ei = E2 in the above equation. 

Magnetic Dipole Interaction 

For the magnetic dipole case, with g = fi and F^,^ = cr^j, in eqn. (1B6[) . the averaged 
amplitude squared reads 



4g4 



{Pi ■P2)[iki ■pi)ih ■P2) + (^2 -Pilih-Pi) + {Pi ■P2)mim2] 



and this equation in the CM frame looks like 

mim2 



[s — m? — 



2E1E2 



E1E2 J 

The total cross section is then given by 

mim2\ 1 A {s — rnl — m^) 



aVrel 



57r 



1 + 



E1E2 J 3 



2E1E2 



(Bll) 



(B12) 



(B13) 



One can obtain eqn. (fT3i) by taking mi = m2 and Ei = E2 in the above formula. 



Appendix C: Event Rate for DM-Nucleus Elastic Scat- 
tering through Scalar Interaction 

Here we consider scattering of DM and nucleus through scalar mediated spin-independent 
interaction. The details of this calculation are given in [2]. Here we shall give only the 
results of the calculation. The differential recoil rate per unit detector mass for DM-proton 
scattering can be written as 



dR p^apA^ 



dE 



R 



dv 



f{v) 



(CI) 
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Figure 8: Plot shows variation of spin- independent DM- nucleus scattering cross section as a 
function of the DM mass for different experiments. 



where q = y/rriNER is the nucleus recoil momentum, ap is the DM-proton cross-section, A is 
the number of nucleons and fip is the reduced mass of the DM-proton system. Here F[Eji) 
is the Woods-Saxon form factor given by 



F(Er) = ^2lMe-rs^ (C2) 
qr 



where q = v^2m7v^E/j is the momentum transferred, m^r being the nucleus mass, s ~ 1 fm 
and r = y/R^ - 5s^ with R = 1.2^^/3 jj^g 

spin-independent differential cross section for 

contact interaction scattering is given by 



Gl-,F\En) = -T^F'iEn) (C3) 



where 



C = —^[Zfp+iA-Z)UY (C4) 



Jo "|qr 

Here |qp = 2mNER is the magnitude square of the momentum tranferred, //at is the 
reduced mass of the DM-nucleus system, Z is the number of protons and A is the number 



20 



of nucleons. Taking fp ~ we can write the DM-nucleus cross section in terms of the 
DM-proton cross section [7] using eqn. OCSp as 

ao = (C6) 



where /ip is the reduced mass of the proton-DM system. Substituting from eqn. flC6p in 
eqn. flCSp we have 

The differential rate per unit detector mass can be written as 

where the quantity / is as defined in the main section. Substituting from eqn. (IC7p in the 
above equation we get 

dR = P^''p^\ f\Er) ( ^) dvd\q\' (C9) 
/I™ ™ .^i y V J 

|2 



Using (i|qp = 2mN dEpi we get the differential rate per unit recoil energy, 

dR _ p^apA^ ^2ir. ^ f /(^) 



Jv>Vmin 



-dv 



dER 2mx/i2 
which is same as eqn. dCip . 
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